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Abstract: This paper is concerned with a novel and efficient algorithm, the pre-sweeping algorithm, 
for solving the Chan-Vese model in image segmentation. This algorithm avoids solving 
the Euler-Lagrange equation corresponding to the energy function. It improves the com- 
putational speed dramatically when applied to the Chan-Vese image segmentation model 
and keeps all the advantages of the level set method such as automatically handling the 
topological variety, preserving the sharp-angle of the curve, etc. Our algorithm is also 
efficient for some image that the Song’s method does not work on when using the Jacobi 
iteration. Moreover, our algorithm is easy to extend to arbitrary finite dimensional image 
segmentation. 
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1 = Introduction 


Image segmentation is a bridge between image processing and image analysis. The aim is 
to partition a given image into its constituent parts with same or similar properties, such as 
intensity, colors, and textures etc. It has been widely used in many fields, for example, the edge 
detection in computer vision, medical image processing, character recognition, and classification 
of remote sensing image!?-5]. 

Variational method is a powerful tool in image processing and other fields such as shape 
optimization, flow control, material science!*!2). The basic idea of the variational method is 
to minimize an energy/cost function depending on the problem to be dealt with. The Kass- 
Witkin-Terzopoulos model!!3! may be the first effort in this direction. Since the energy function 
in the Kass-Witkin-Terzopoulos model is nonconvex, no uniqueness result is available and the 
result is much depending on the original contour. Mumford and Shah presented a model!"4] with 
good theoretical properties. The model has difficulty to be applied since it requires to minimize 
the energy function with two variables in two different natures. Moreover, these models identify 
boundary with the edge detector. So for the edges which is not defined with gradient, these 
methods are poor efficient. 

To improve the efficiency of the Mumford-Shah model, Chan and Vese!!! proposed a new 
method introducing level set!9:15.16) into the Mumford-Shah function!"4] for image segmentation. 
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Experience has shown the efficiency of C-V model for image segmentation. 

Many new and efficient methods have been proposed to solve the Chan-Vese model in recent 
yearsl!7-20]_ Among these methods, the fast algorithm proposed by Song and Chan!"4] is a novel 
and efficient one. And it can be applied to broader range of optimization problems without 
solving PDE. We propose an algorithm which is similar but different from the Song’s algorithm. 

All the examples considered in this paper are 2-D gray images. It is not meant the method 
is limited to 2-D images. It is convenient to extend our algorithm to arbitrary finite dimension. 

This paper is organized as follows. In the next section, we briefly review the C-V model. In 
section 3, we introduce a new algorithm, pre-sweeping algorithm, for C-V model. In section 4, 
we give several examples of our algorithm applied to image segmentation. And we compare it 
with the algorithm of active contour method (ACM) and the fast algorithm proposed by Song 
and Chan!"8], Section 5 is the conclusion and some remarks. 


2 Chan-Vese model for image segmentation 


The C-V model is a variational model for 2-phase image segmentation. From now on, let 
Q C R? be a bounded open set. Curve C(s) = (z(s), y(s)) : R — R? divides Q into two regions: 
Dı (inside of C(s)) and Dz (outside of C(s)). Dı and Dz are open sets. up : R? — R is 
an original image. A piecewise constant function u(x) is a segmentation result of up. u(x) = 
ei(z € Dı) and u(x) = ce(x € D2). Then the energy function is 


Fov(C,¢i,¢2) = a- (length(C))” + 8- area( Dı) 
+À: | Juo — cı |?da + ` f Juo — c2|°dz, (1) 
Dı D2 


where a, 8, Àı and Àz are weight parameters. The parameter can be given in advanced accord- 
ing to the image to be processed. In the C-V model, cı and cz are mean values of uo(z) on the 
region of Dı and Dg, respectively 
u(x)dx u(z)da 
Cı = Jp wields ; C2 = Jpg te )de X (2) 
meas( Dı) meas( D2) 
Introducing the Heaviside function H(z) = 0(z < 0), H(z) = 1 (z > 0) and a Lipschitz 
continuous level set function ¢(z,t) : R? x [0,co] > R 


>0, «ED, 
d(z,0)¢ =0, rE, (3) 
<0, xe Dz, 


into the energy function. We use the zero level set of the level set function ¢(z) : R? — R to 
represent the curve C(x). So the energy function can be written as|!+19,20] 


F(ġ,c1,c2) = a( f | H(¢) jae)" +8 f Hae 


+21 i luo — c1 |H (ġ)dz + af luo — c2|° (1 — H(¢))dz. (4) 
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The classical way to solve the problem 
ae F (4, cj ($), c2(¢)) 


needs to solve the corresponding Euler-Lagrange equation, usually a nonlinear parabolic PDE. 
Due to the Courant-Friendrichs-Lewy (CFL) condition, the numerical algorithm is often not 
efficient. Song and Chan!'8) proposed a fast algorithm without solving PDE. The algorithm is 
efficient and converges very fast. Motivated by this algorithm, we propose an algorithm that is 
similar but different from the Song and Chan’s algorithm. 


3 Pre-sweeping algorithm for Chan-Vese model 


We suppose that the level set function ¢(x) = 0 divides Q into two different regions: D, 
and D2. Let kı and kg be the number of pixels in Dı and Dg, respectively, and cı and c2 be 


mean value of u(x) in the Dı and D2. Let 
d(x) =1(x € Dı), G(x) =-l(x€ D2), Do = {| G(x) =0; 2 € O}. 


We use ko and cg to present the pixels number belonging to Dp and the mean value of Do, 


respectively. 
Let H(z) =1(z> 0) and H(z) =0(z <0). We rewrite the energy function the Chan-Vese 


model as follows 


F(¢(z2)) = a f | uo — cı KONEEN | uo — C2 |? (H(—¢))dz 
i Í E T ES) T | VH(d) ldz. (5) 


Where f Do fdz represents the integral on Do. Do may be a curve or a surface. Then the image 
segmentation problem is to minimize above function in piecewise constant function (PCF) 


min E (9). (6) 


Set the parameters Ag = Ay = A2 = 1 and u = 0. Let AF, and AF; be the varieties of 
the energy function (5) when we change a point 29 € Do to Dı and to D3, respectively. If we 
change ġ(zo) = 0 to ¢(zo) = —1, the new energy is 


Pl) = X ola) -a+ Y fuolz)-&P+ J hoe) -êl (7) 





zEeD: xzE(Dz2U zo) zE(Do\zo) 
where 
7 k2c2 + uo(£o) NE koco — uo{xo) 
nm feet OO ko— 1 


So the difference between the new energy and the old energy is 


- (ule) 0)" 5.8) 


AF (z0) = F2(x0) F(z) = (uo(xo) = c2) 








2 k2 
ko +1 
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Similarly, it is easy to calculate 








AF (z0) = (uo(zo) = a) I = (uo(z0) = oo) 1 (9) 


Our pre-sweeping algorithm for solving the problem (6) is: 

Initialize For a given image ug and an initial curve ¢(x), we calculate the cı and c2. 
Without loss of generality, we assume cı > c2. Set d(x) = 1 if z € Di, d(x) = -1 if z € D2, 
and Do = 2: (x) = 0. 

Step 1 Pre-sweeping. For every x € Q, if d(x) = 1 and u(x) < c2, then change ¢(x) = —1; 
else if c2 < u(x) < cı, change (x) = 0. If (x) = —1 and u(x) > cı, then change ¢(x) = 1; 
else if c2 < u(x) < c1, change (x) = 0. 

Step 2 Sweeping. For every point x € Do, compute AF; (x) and AF (x) with (8) and 
(9). If AF, > 0 and AF» > 0, then stop; otherwise change ¢(x) to 1 or -1 according to the less 
AF;(xz). Then go to the beginning of Step 2. 

Theorem 3.1 Ifup:Q— R is a 2-phase image and satisfies 

a, TEA, 
uo (x) = (10) 
b, zEB, 
where a and b are constants (a > b). Then for any initial curve ọ(x) = 0 dividing Q into two 
different regions, the pre-sweeping algorithm converges to the exact result after pre-sweeping 
for the model (5) with ào = A1 = à2 = 1 and p = 0. 

Proof Suppose that the curve (x) = 0 divides Q into Dı and D2. cı and cz are the mean 
values of Dı and De, respectively. Without lose of generality, we assume cı > c2, then we have 
b< cz <cı <a. Set o(z) = 1 if z € Dı and o(x) = —1 if z € Do. 

If z € Dı and uo(£) < c2, so we have uo(xz) = b. According to the algorithm, we change 
(x) = 1 to d(x) = —1. If x € Dz and uo(z) > ci, so we have uo(x) = a. Using the pre-sweeping 
algorithm, we change ¢(x) = —1 to g(x) = 1. 

So after the pre-sweeping, (x) = —1 for x: u(x) = b and (x) = 1 for x: uo(z) = a. 

After pre-sweeping, we sweep on the Do by using the following theorem: 

Theorem 3.2 Suppose Do # Ø. Let 


A= {a | u(zo) < u(r) < 2co — u(zo)}, B = {z | u(x) > u(x) > 2co — u(zo)}. 


1) If zo is a solution of mi {u(x)}, and the energy function (5) is non-decreasing when 
changing $(xo) = 0 to ¢(ap) = —1, then the energy function is non-decreasing for every z € A 
changing to ¢(x) = —1. 

2) If Zo is a solution of max{u(z)}, and the energy function (5) is non-decreasing when 
changing (zo) = 0 to ġ(zo) = 1, then the energy function is non-decreasing for every x € B 
changing to ¢(x) = 1. 

Proof We only need to prove the 1). 2) can be proved in the same way. 

1) The energy function (8) is non-decreasing when changing (zo) = 0 to (zo) = —1, so 
AF (xo) > 0. For every x € A, we have 


c2 < ulto) < co, [eo — u(zo)] > u(x) — col, lez — u(z0)| < |u(z) — cal. 
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So we have 
k l> 
AF(a) = (u(z) c2)” i (u(x) co)? E i 
k 
> (u(x) J a) i = (u(xo) = co) 1 > 0. (11) 


In the above theorem, we have not considered the x : u(x) > 2cg — u(zo) when changing 
$(zo) = 0 to $(xo) = —1 and the x: u(x) < 2c9 — u(xo) when changing (a9) = 0 to ġ(zo) = 1. 
It is obvious that we want to move the x with lager value to the domain Dy, and the small value 
to D2. So the method is reasonable. 


4 Numerical experimental result 


In this section, we present several examples to show the application of our pre-sweeping 
algorithm for image segmentation. For simplicity, we set Ay = A; = Àz = 1 and & = 0 in all 
the examples. The pre-sweeping algorithm is very easy to program. All of the examples in 
the paper are executed in a personal computer with CPU Pentium IV 3.40G Hz and 512M 
memory. We program with Matlab. 

In Figure 1 and Figure 2, we applied three different algorithms, our pre-sweeping algorithm, 
the Song and Chan’s fast algorithm!"®), and the active contour method without edge!"), to an 
image named cameramen with the same original level sets all being given as a circle. This 
example shows that our algorithm and Song’s algorithm are almost convergence to the same 
result. On the other hand, we will show the proposed algorithm is more efficient than the fast 
algorithm proposed by Song’s later on. 

In Figure 2, we use the active contour without edge method to the same image. The 
result shows that the active contour without edge method converges very slowly. We show the 
segmentation results iteration in 1000, 6000 and 10000 iterations. And it requires almost ten 
thousands times iteration if one wants to get the similar result as the former two methods got. 





(a) (b) (c) (a) 


Figure 1: (a) The cameramen image with the initial level set; (b) and (c) are segmentation 
results of the cameramen image using the pre-sweeping algorithm and the Song’s algorithm, respec- 


tively; (a) The result of level set using pre-sweeping algorithm. 
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(a) (b) (c) (a) 


Figure 2: (a) The image with the initial level set; (b), (c) and (d) are the level sets after 1000, 
6000, and 10000 iterations, respectively. 


Table 1 compares our pre-sweeping algorithm with the Song’s algorithm in terms of CPU 
times. We apply the two algorithms to the same image with different sizes. The bigger size the 


image is, the more advantage our algorithm shows. 


Table 1: Comparison of CPU times 











image size 128x 128 256x256 512x512 1024x1024 
pre-sweeping algorithm 0.1879s 0.2399s 0.5278s 1.6844s 
Song and Chan’s algorithm 0.2506s 2.5108s 30.8046s 237.3616s 





In Figure 3, we give an image which the Song’s algorithm with the Jacobi iteration does 
not work. If the initial level set function is 6 = 1 on the left side and ¢ = —1 on the right 
side, using the Song’s method with Jacobi iteration, when a point z in left can be changed to 
(x) = —1 or d(x) = 1, the point on the right corresponding would also be changed to ¢(x) = 1 
or ¢(x) = —1 for symmetry. So the Song’s algorithm failed for this kind of images. While our 
pre-sweeping algorithm can converge to the accurate result whatever using the Jacobi iteration 


or the Gauss iteration. 





(a) (b) 


Figure 3: (a) A type of symmetrical images that the Song’s algorithm does not work when using 


the Jacobi iteration; (b) The segmentation result of our pre-sweeping algorithm. 
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In figure 4, we use our algorithm to segment an image with three different initial conditions. 
All of them converge to the accurate result. This shows our algorithm is quite robust. And it 
can be seen that the interior contour of the objects can be detected automatically. Even for 


the non-convex contour, it can also be detected correctly. 





(a) (b) (c) (d) 


Figure 4: To segment an image using our pre-sweeping algorithm with three different original 


conditions (a), (b) and (c); all of them get the same result which is shown in (d). 


In Figure 5, we show an image with three different objects. It requires two level sets to 
divide the image correctly when using the Song’s algorithm. Our algorithm only needs one 


level set to get the correct result. 





(a) (b) 


Figure 5: We use one level set to divide an image into three different parts. (a) is the original 


image with initial level set; (b) is the result. In this case, the level set @ = 0 represents the third part. 


The results of our algorithm used to the noised image are shown in the F igure 6. We use the 
pre-sweeping algorithm to the image with 5 percent salt and pepper noise. The two results are 
shown together with the same initial level set. From the figures we can see that our algorithm 


can get satisfied segmentation result even for the image with noise. 


5 Conclusion and some remarks 


We proposed a new algorithm, the pre-sweeping algorithm, in this paper. The new algorithm 
is very efficient and robust. It has many advantages compared with other algorithms used to 
solve the Chan-Vese model in image segmentation. The algorithm preserves the advantage of 
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the level set method. Our algorithm doesn’t require to solve any PDE. So using our algorithm, 
it doesn’t need to consider the CFL condition. This makes our algorithm to converge faster 
than the PDE based methods. For some 2-phase images, our algorithm only needs pre-sweeping 
to get satisfied result. And for some images which have three different parts, the pre-sweeping 
algorithm only needs one level set to get correct result. 

All the examples in our paper are 2-D images, but the algorithm can be extended to high 
dimension image segmentation easily. Moreover, the algorithm can be used to a lot of opti- 
mization problems when the energy function could be easily calculated. 

Since we need not calculate the total energy, the length term in (5) is ignored in our algo- 
rithm. So the method is sensitively for noised image. But we can get satisfactory result using 


the pre-sweeping to the denoised image. 





(a) (b) (c) 


Figure 6: (a) The original image added with 5% salt-and-pepper noise with initial level set; (b) The 
segmentation result of pre-sweeping algorithm to the noisy image directly; (c) The segmentation result 


of the image after average filtering. 
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